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Abstract 

We study the decays i% — ► i\HH, — > i\ZZ, I2 — > t\W + W~ and £2 — > t\ZH, with an 
eye towards measuring the CP-violating supersymmetric parameters contained in these 
processes. We find that I2 — » t\HH tends to have the largest CP asymmetry and width, 
and is perhaps the most favourable experimentally. These decays are sensitive primarily 
to 4>A t > the phase of the trilinear coupling A t . 
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Supersymmetry (SUSY) is widely thought to be the physics that lies beyond the 
standard model (SM). SUSY theories typically contain many new parameters, some 
of which are complex, and hence violate CP. Should SUSY be found experimentally, 
we will want to find the values of all of these parameters. In particular, it will be 
important to measure the CP-violating SUSY phases. 

CP- violating SUSY effects have been studied extensively in meson mixing [1], 
in CP violation related to the 5-meson system [2] and in electric dipole moments 
(EDMs) [3]. In particular, EDMs provide quite stringent constraints on the low- 
energy CP-violating SUSY phases of the superparticle couplings. For example, if 
sfermion masses are taken to be of order the weak scale and complex SUSY parame- 
ters have phases of order unity, theoretical predictions for EDMs are not generically 
in agreement with experimental limits. Nevertheless, the so-called SUSY CP prob- 
lem can be avoided in several SUSY scenarios [4] as calculations of EDMs are highly 
model-dependent. In this paper, all SUSY parameters used as inputs in the numeri- 
cal simulations are assumed not to violate EDM constraints. The approach adopted 
here is that our observables may offer an independent measurement of relevant CP- 
odd SUSY parameters. 

In a recent paper [5], we showed that the decay t 2 — > tir~T + is particularly 
sensitive to 0^ t , the phase of the trilinear coupling A t . (The "stops" i\ and t 2 are the 
two mass eigenstates of the scalar superpartners of the top quark, with mf 2 > m^.) 
In the present paper, we examine how some other decay processes depend on the 
SUSY phases. Throughout we assume that SUSY has been discovered, and that the 
CP-conserving parameters (e.g. masses of SUSY particles) are known independently. 

All effects which violate CP require the interference of (at least) two amplitudes. 
For a given decay, there are two types of CP-violating signals. Direct CP asymme- 
tries (spin-independent or spin-dependent) are proportional to sin S, where S is the 
relative strong (CP-even) phase between the interfering amplitudes. Triple-product 
(TP) asymmetries take the form V\ ■ (v 2 x ^3) (each Vi is a spin or momentum), 
and are proportional to cos 5. One can therefore have a nonzero TP even if S — 0. 
It is also possible to have a nonzero TP with only a single decay amplitude if, for 
example, the intermediate particle has both scalar and pseudoscalar couplings [5, 6]. 

Strong phases can be generated in one of two ways. First, one can have the 
exchange of gluons between the particles involved in the decay, leading to QCD- 
based strong phases. Unfortunately, we do not know how to calculate the strong 
phases in this case. Alternatively, the strong phases can be generated by the (known) 
widths of the intermediate particles in the decay. Given that we want to measure 
the SUSY CP phases, and not simply detect the presence of CP violation, we must 
consider decays in which the strong phases are generated in the second way. Thus, 
the decay processes cannot contain too many particles which couple to gluons. 

It is therefore quite natural to consider the decays t 2 — > t\HH, t 2 — > t\ZZ, 
t 2 — > t{W + W~ and t 2 — > t\ZH, where if is a neutral Higgs boson. There are several 
points here which we should note in relation to these processes. First, SUSY theories 
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involve two Higgs doublets which contain (in the gauge basis) two neutral scalars 
and one pseudoscalar. In the mass basis, these particles mix and one obtains three 
mass eigenstates Hi, H 2 and H 3 . All four decays receive contributions from several 
diagrams, including those with an intermediate Hi (i = 1,2,3). The amplitudes 
are typically dominated by those diagrams in which an intermediate particle can 
go on shell. In SUSY theories, the lightest mass is m Hl = 0(100) GeV and we 
therefore take the final-state H to be H x . When it appears as an internal line in 
a diagram, the Hi is too light to decay (on shell) to H^Hi, ZZ, W + W~ or ZH 1 . 
Nevertheless, we shall retain such diagrams since they can in principle give non- 
negligible contributions. 

Second, triple products are due to terms of the form Tr [7^7^^70-75] in the square 
of the amplitude. Since no fermions are involved in the decay processes, no CP- 
violating TP's can arise here, and we have only direct CP asymmetries. The final- 
state particles do not couple to gluons and any exchange of gluons between t 2 and t\ 
only serves to renormalize the couplings of the stops. Thus, the strong phase arises 
only due to the widths of the intermediate particles H 2 and H 3 . For a given set 
of SUSY parameters, the widths r 2 and r 3 are calculable (as are the 'off-diagonal 
widths' associated with transitions Hi <-> Hj [7]). Thus, the measurement of CP 
violation in these decays (due to direct CP asymmetries) will allow us to extract 
and/or constrain the SUSY parameters, including the CP- violating SUSY phases. 

Third, there can be significant interference between the two decay amplitudes 
only if the masses of H 2 and H$ are similar. Fortunately, it is relatively common in 
SUSY theories that m# 2 ~ m# 3 [8]. Since it is assumed that the masses are known, 
we will know beforehand whether or not CP violation is likely in these decays. 

Finally, we make a comment regarding our previous work, which examined CP 
asymmetries in i 2 — > i\T~T + [5]. The process i 2 — > tiT~T + is interesting since it is 
theoretically clean (no strong phases from gluons) and it can have large CP asymme- 
tries while simultaneously having a moderately large branching fraction. This pro- 
cess is particularly attractive if the intermediate Higgs bosons are too light to decay 
to heavier final states, such as those considered in the present work (PU + 1U~, H 1 H 1: 
etc.). Nevertheless, as noted in Ref. [5], the simple rate asymmetry for t 2 — > t\T~r + 
is extremely small. Thus, sizeable CP asymmetries are only expected in cases where 
the spin of one or both of the final-state leptons is measured. The processes con- 
sidered in the present work have an advantage over t 2 — > i\T~T + in that they do 
not require the measurement of any spins - the regular rate asymmetries can have 
relatively large values. 

There are five classes of diagrams that contribute to the four processes under 
consideration, although not every class of diagram contributes to each process. The 
first class is shown in Fig. 1. In this case the heavier stop decays to the lighter 
stop and emits a Higgs boson. The Higgs then decays to the final state f\f 2 (where 
/1/2 = H\H\, etc.) This process can proceed resonantly if the intermediate Higgs 
boson(s) can go on shell. Figure 2 shows the other four possibilities. In diagram 
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FIG. 1: Feynman diagram for the decay t 2 — ► t 1 /1/2 with /i/ 2 = H\Hi, ZZ, W + W and 
Zi^. The Higgs propagator has off-diagonal terms, so transitions Hi — > iJj are allowed. 



(a) the stop decays to particle / 2 and a squark ^ (either a stop or a sbottom), 
which subsequently decays into j\ and the lighter stop. Diagram (b) is a crossed 
version of diagram (a). Diagram (c) is similar to that shown in Fig. 1, but with the 
intermediate Higgs replaced by a Z. Finally, diagram (d) shows the contribution 
due to a possible quadrilinear vertex. 

Of all the diagrams shown in Fig. 2, there is only one case in which the process 
can proceed resonantly. This occurs for the decay t 2 — > iiW + W~ (/1 = W + and 
f 2 = W~), in which case the intermediate squark in diagram (a) is a sbottom. 
When considering t 2 — > t 1 W + W^ in our numerical work we choose parameters in 
such a way that the sbottom(s) cannot go on-shell. Thus, in all cases the diagrams 
in Fig. 2 contribute non-resonantly. While one might be tempted to ignore these 
diagrams compared to the resonant Higgs contributions in Fig. 1, we have found 
that the non-resonant diagrams can sometimes give non-negligible contributions, so 
we include them in our calculation. We will discuss this point further below. 

We turn now to a calculation of the width and rate asymmetry for the process 
^2 —> ^1/1/2, with fif 2 = H1H1, ZZ,W + W~ and ZH\. We define two invariant 
masses as follows: 

M 2 = (p 1+ p 2 ) 2 , (1) 
P 2 = (Pi+Pt,) 2 , (2) 

where p± t 2 are the four-momenta of f\ and / 2 , respectively, and is the four- 
momentum of the t\. All dot products that arise in the calculation may be written 
in terms of M 2 , p 2 and the various particle masses. 

As noted above, there are no TP's in this process. In principle there are 
polarization-dependent CP-violating observables similar to the single-spin CP asym- 
metry defined in Ref. [5]. It turns out that, in the limit in which one can neglect 
the non-Higgs-mediated diagrams, such polarization-dependent observables are sen- 
sitive to the same combinations of underlying SUSY parameters as are the rate 
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(c) (d) 

FIG. 2: Feynman diagrams that can contribute to the decays — > ff/i/b with /1/2 = 
H\Hi, ZZ,W + W~ and ZH\. Not all diagrams contribute to each process. 



asymmetries. Thus, in this limit, one gains no new information by measuring po- 
larizations. In principle the cross terms of resonant and non-resonant amplitudes 
can give new contributions to polarization-dependent CP-odd observables, but such 
contributions are expected to be suppressed. Given the difficulty of the related mea- 
surements, and the possible suppression, we will ignore such observables and sum 
over the polarizations of the final-state particles. 

Let us consider first the "processes" — > ff/i/2 ( a s opposed to the "anti- 
processes," involving the decay of a t%, which will be considered in a moment). (Note 
that the indices ± associated with the i's indicate that they have charge ±2/3.) The 
Higgs- exchange diagrams shown in Fig. 1 contribute for all four final states. The 
following are some technical details for each process regarding the diagrams shown 
in Fig. 2. 

1. — > iiHiHi. Diagrams (a), (b) and (d) in Fig. 2 contribute. (The H^-H^-Z 
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coupling is zero.) The intermediate <jj are stops. These stops cannot go on 
shell. 

2. — > t\ZZ: Diagrams (a), (b) and (d) contribute. The intermediate q~j in 
diagrams (a) and (b) are stops; they cannot go on shell. 

3. — > tiW + W~: Diagrams (a), (c) and (d) contribute (fi = W + and f% = 
W~). The intermediate squarks in diagram (a) are sbottoms. In principle the 
sbottoms could go on shell, although in practice we choose parameters such 
that this does not occur 1 . 

4. — ¥ t\ZH\. Diagrams (a), (b) and (c) contribute. The intermediate top 
squarks in diagrams (a) and (b) cannot go on shell. 

The amplitudes for the four processes may be written in the following manner: 

Ahh = B H h , (3) 

Azz = (BzzgT + CzzV^vl + D zz^ tl V\ + EzzP&i) e^ 2 * , (4) 

Aww = {Bwwgv + C wwP £pl + D wwV ^ V \ + e^ 2 * , (5) 

Azh = (Bzh^ + CzhpQ^, (6) 

where the e^i?* are the polarization tensors for the final-state vector mesons. The 
parameters B^f 2 , . . . , E^f 2 are functions of the various coupling constants and mix- 
ing matrices. In each case, B^f 2 contains the Higgs-mediated contributions, as well 
as one or more other contributions coming from diagram(s) in Fig. 2; that is, 

B hh = Bj;r + 6B flh , (7) 

where 5Bf l f 2 denotes the non-Higgs-mediated contributions. Expressions for 5Bf x f 2 , 
C/1/2) -^/i/2 an< i Efih ma y ^ e f° un( i m the Appendix. 

We include here the expressions for the Higgs-mediated contributions, Bj^ s , 
since they are of the most interest to us (they tend to dominate, and they are also 



If one or both sbottoms go on shell the calculation becomes more complicated. In particular, 
the asymmetry would depend on the total widths of the sbottoms. Furthermore, gluon exchange 
between ti and t\ in diagram (a) could give rise to additional strong phases. Even if one were to 
include the on-shell sbottoms diagrams, the resulting asymmetry would tend to decrease. The 
sbottom diagrams are unlikely to interfere with each other, since sbottom masses are generically 
not close to each other. Also, since sbottoms do not proceed via the s-channel, the on-shell 
interference between Higgs bosons and sbottoms is restricted to a very small region in the Dalitz 
plot. 
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the source of the required strong phases). These contributions are given by, 

B%T = -fY,9H&hDii(M*)9n i H l Hjl i , (8) 



h3 



D Higgs vgm w x - (t\j2\ / n \ 

B zz = ~ C0S 2 dw l^dH^kDijiM )g HjVV , (9) 



B ww = -v9m w ^2g H . { * il Di j (M 2 )g H . vv , (10) 



5 ^ SS = ^E^&^W^z , (ii) 

costv ^. 

where we have included factors rji = 6 and 772,3 = 2 in order to correctly account 
for the manner in which the gHiHjH k are defined in Ref. [9]. The Higgs propagator 
matrix is given to a good approximation by [7, 10] 

iD(M 2 ) = 

I M 2 — m 2 Hi + i Imflii i Imft^ i Imfli3 \ 

ilmfl 21 M 2 - m 2 H2 + ilmfl 22 ilmU 23 ^ (12) 

V i Imfiai i Imn 32 M 2 - m 2 H3 + i lmf[ 33 J 
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Expressions for the absorptive parts of the Higgs-boson self-energies, Imllj^M 2 ), 
may be found in Ref. [7]. The Appendix of the present work contains a brief discus- 
sion of the various couplings g Hi i*^ g^mm, Qh.vv and g Hl H 3 z- Of these couplings, 
only gjj.i*^ is complex. We adopt the notation of Ref. [9] for these couplings, except 
in the case of g Hl H 3 z- 

To obtain the width for t 2 — > t± f\f 2 , we multiply the respective amplitude by 
its complex conjugate, sum over the polarization states of the vector boson(s) in 
the final state (if appropriate) and integrate over the squares of the two invariant 
masses, M 2 and p 2 , to obtain 

r «' = 256l^/(?l^| 2 ) W ' < 13 > 

tl \polns / 

where S F = 1 for fj 2 = WW,ZH 1 and S F = \ for fj 2 = E X H X ,ZZ. For 
the H\H\ case there is no sum over polarizations and the calculation is relatively 
straightforward. For the ZZ and WW cases there are two sums over polarization 
states (one for each vector particle in the final state), leading to ten separate terms, 
each with its own kinematical factor. The terms are proportional to |£?| 2 , Re (BC*), 
etc. The ZHi final state only requires one sum over polarization states, resulting in 
three separate terms. The appropriate range of integration for M 2 and p 2 may be 
found, for example, in Ref. [11]. 
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The width T f x f 2 can be similarly defined for the CP-conjugate process t\ — > 
tffif2- 2 The only difference compared to the expression in Eq. (13) is that one must 
complex conjugate the weak phases in B^^, . . . ,Ef 1 f 2 , which amounts to making 
the replacements gn t t*t k ~^ 9*h I*1 anc ^ ^* ^ * n ^ ne var i° us expressions, as well 

as « — * — i in Eq. (11). Some further discussion may be found in the Appendix. ([/* 
is the stop mixing matrix.) For the rate asymmetries to be non-zero, we require (at 
least) two interfering amplitudes, and these must have a non-zero relative strong 
phase as well as a non-zero relative weak phase. The weak phases appear in the 
various couplings and the strong phases are provided by absorptive pieces in the 
Higgs propagator matrix. (In principle, there are also strong phases associated with 
the widths of the Z and the squarks in the diagrams in Fig. 2. The effects of these 
widths are suppressed, however, since the associated (s)particles are always off-shell 
in our calculation.) 

The rate asymmetries for the four cases are then defined to be 

A CP (t2 - fi/i/ 2 ) = ^ /l/2 7^ /l/2 • (14) 

Our previous paper [5] contains an extended discussion of the behaviour of CP 
asymmetries as functions of the invariant mass M. The analysis in the present case 
is complicated by the non-negligible contributions of diagrams that are not in the 
s-channel. Nevertheless, one could still define differential widths as functions of M, 
and these would still be expected to exhibit resonant peaks near M ~ fnu^- 

We now turn to a numerical investigation of rate asymmetries for t 2 — > 
tiHiHi, tiZZ, tiW + W~ and tiZH 1 . In our numerical work we have made extensive 
use of the computer program CPsuperH [9]. Note that since we are using the widths 
of the Higgs bosons to provide the required strong phase, the rate asymmetries will 
only be non-negligible if there is a significant overlap of the interfering resonances 
(see Ref. [5] for further discussion of this point). Fortunately, as noted above, it is 
not uncommon to have m# 2 fa in SUSY. 

Figure 3 shows several scatter plots of the rate asymmetries for t 2 — > 
tiHiHi, t\ZZ, tiW + W~ and t\ZH\. In these plots we have allowed several SUSY pa- 
rameters to vary in specified ranges, taking m#± G (160 — 500) GeV, tan/3 G (1 — 15), 
fi G (200 — 1400) GeV (without loss of generality, the \x parameter is taken to be real 
and positive - see Ref. [5] for further discussion), mg 3 , m^ 3 , m^ 3 G (300 — 700) GeV, 
| A t | G (400 - 2000) GeV and <f> At G (0° - 360°). All other input parameters have 



2 Since the Higgs bosons are typically mixtures of scalar and pseudoscalar states when CP is 
broken, some of the final states that we consider are not eigenstates of CP. Nevertheless, the 
widths we calculate for what we are calling the "CP-conjugate" processes are in fact the widths 
that would be measured experimentally. Furthermore, the asymmetries that we calculate are 
zero in the CP-even limit. 
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FIG. 3: Scatter plots showing rate asymmetries for I2 — > i\H\Hi, t\ZZ, t\W + W~ and 
t\ZH\. The horizontal axis in each plot gives the sum of the widths for the process and 
the anti-process for the decay in question. 



been assigned fixed values and have been taken to be real. 3 We have insisted that 
the mass of the lightest Higgs boson be greater than or equal to 100 GeV and have 
chosen values for the charged Higgs mass and tan (3 that are consistent with the 
recent bound from Belle, tan(3/m H ± < 0.146 GeV -1 [12]. We have also insisted 
that all supersymmetric particles have masses greater than or equal to 100 GeV 
and that the lighter stop be kinematically allowed to decay by at least one of the 
following two modes: t\ — > tXi or tf — > bxf- For the purpose of these plots we have 
set the widths of intermediate stops and sbottoms to be 10 GeV [13]. (Varying these 
widths does not seem to have a significant effect on the numerical values obtained 



3 For completeness, we list here some of the other parameter choices, in CPsuperH notation: 
Mi = 100 GeV, M 2 = 200 GeV, = M 3 = 1000 GeV, m L ^ = 150 GeV, = 600 GeV, 
Ab = 1000 GeV and A T = 750 GeV. With these parameter choices, and after applying cuts noted 
subsequently in the text, some of the supersymmetric particles have the following mass ranges: 
mo T G (136-138) GeV, m h g (199-558) GeV, m t ~ 2 G (468-883) GeV, m Si G (255-685) GeV, 
m fl G (146 - 156) GeV, m~± G (196 - 200) GeV and m^o g (100 - 101) GeV. Also, the mass 
ranges for the three Higgs bosons arc G (100 — 126) GeV, itih 2 g (161 — 490) GeV and 

m H3 e (161 - 495) GeV. 



8 



for the asymmetries, as long as the sbottoms are not allowed to go on shell.) The 
horizontal axes in these plots show the sum of the widths for the process and the 
anti-process. If we assume a typical stop width to be approximately 10 GeV, we 
see that it is possible to have large rate asymmetries while simultaneously having 
relatively large branching ratios. For example, in t 2 — > t\H\H\ it is possible to have 
an asymmetry with magnitude of order 20—30% when T HlHl + T^Hx — 0.2 GeV. 
Among the decays considered, i 2 — > tiH\Hi is clearly favoured in that it tends to 
have a larger width for a given CP asymmetry. Experimentally, however, some of 
the other decay channels may be favoured due to increased detection efficiencies. 
Note that the WW final state has fewer data points compared to the other plots 
because parameter sets for which a sbottom would have gone on shell have been 
discarded for that plot. 

In Figure 3 we have allowed the phase of the trilinear coupling A t to vary, but 
have set all other phases to zero. In addition to the results shown in Fig. 3, we have 
also performed a limited analysis of the effects on the CP asymmetries due to the 
phases of other SUSY parameters. The strongest effects come from <pA b and 0m 3 , the 
phases of the trilinear coupling Ab and gluino mass M 3 , respectively. The phase 4>A b 
can affect the asymmetries due to the involvement of A^ in the mixing of the Higgs 
bosons. The phase of M 3 comes into play through its effect on the Higgs-bottom- 
bottom effective vertex 4 . The phase 4>A b only seems to produce a non-negligible 
effect when \A b \ is greater than 5 TeV (due to the difference in the top and bottom 
Yukawa couplings). Allowing <pM 3 to assume a nonzero value can lead to changes in 
the CP asymmetries of order ±0.1 if IM3I is taken to be approximately 1 TeV. 

To summarize, many physicists believe that supersymmetry (SUSY) is the physics 
that lies beyond the standard model, and that it will be found at a future high-energy 
collider. Assuming that this is the case, we will want to measure the CP-violating 
SUSY parameters. With this in mind, in a previous paper [5], we studied the decay 
i 2 — > tiT~r + . We found that the CP asymmetries could be large, but that they 
require the measurement of the spin of one or both of the final-state r's. 

In the present paper we have studied the decays i 2 — > tiHH, t 2 — > tiZZ, t 2 — > 
t{W + W~ and t 2 — > t\ZH , constructing rate asymmetries that are sensitive to CP 
violation in the underlying SUSY theory. There are several ways in which the current 
work complements that performed in our previous study. First, as we have shown, 
the CP asymmetries considered here do not require the measurement of any spins - 
the rate asymmetries alone are measurable and they can be relatively large. Second, 
we had found in our previous work that t 2 — > t 1 T~T + is sensitive primarily to (f>A t , 



4 Even though the decays under consideration do not involve this effective vertex in a direct way, 
the strength of this vertex can affect the widths of the two heavier Higgs bosons, whose values 
can undergo large variations for certain values of 0m 3 - This could lead, for instance, to a larger 
overlap between the two heavier Higgs bosons in the available phase space, giving rise to an 
enhancement of the CP asymmetries. 
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the phase of the trilinear coupling A t . As noted above, the decays considered in 
this work are also sensitive primarily to <f>A t , although the phases associated with 
the trilinear coupling Aj, and the gluino mass M 3 can also affect the asymmetries. 
Third, it is worth noting that the decays considered here tend to involve a different 
range of Higgs masses than those considered in our previous work. For example, in 
the present case we need to have m^ 23 > 2mw in order for t 2 — > tiW + W~ to yield 
an appreciable asymmetry along with a moderately large branching fraction. In our 
former work, non-negligible asymmetries and branching fractions could be obtained 
for lower Higgs masses. 

Of the processes studied in this work, we find that t% — > t\HH tends to have the 
largest CP asymmetry and width, and so it is the process which is perhaps most 
amenable to study. (Note that experimental detection efficiencies may make other 
decay modes more favourable.) As we have emphasized, the decays of stops depend 
primarily on cf>A t , and so do not provide a sensitivity to a large number of phases. It 
would be of considerable interest to investigate CP asymmetries in analogous decays 
of sbottoms or staus to assess their sensitivity to various SUSY phases (4>A b , 4>a t) 
and probably also <f>A t )- We are currently beginning such a study to determine if 
such decays might provide complementary tools for the study of CP violation in 
SUSY. 
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APPENDIX: EXPRESSIONS FOR COUPLING CONSTANTS AND AM- 
PLITUDES 

In this Appendix we describe some of the Higgs couplings in Eqs. (8)-(ll) in 
more detail. We also provide the full expressions for the amplitudes referred to in 
Eqs. (3)-(6). 

The couplings between the stops and Higgs bosons are defined as follows [9], 

^4 = (r^% 7 O m ^c/i, (A.i) 

where O and U* denote the Higgs and stop mixing matrices [9], respectively, and 
where i = 1,2,3 and j, k = 1,2. As usual, v — Jv\ + v\ is defined in terms of the 



10 



vacuum expectation values of the two Higgs doublets. Expressions for the 2x2 
matrices T at 1 may be found in Appendix B of Ref . [9] . These matrices depend on 
the SUSY parameters A t , /i, cos (3 and sin/3 (where tan/3 = v 2 /v±). It is interesting 
to note that the couplings involving scalar Higgs bosons are real and those involving 
the pseudoscalar Higgs are purely imaginary. These couplings are generally complex 
if CP is broken. 

For the Higgs-Higgs-Z vertex we adopt the same notation as in Ref. [14] (this 
differs somewhat from that employed in Ref. [9]): 

Chhz = o^&a ^ dHiHjzZ^Hi Hj , (A.2) 

Z COSt/ w j>i 

where 

= f 3i (cp0 2j - spOij) - 3j (cp0 2i - spO u ) , j > i 
yHiHjZ j Qj otherwise, 

with cp = cos (3 and sp = sin/3. The Higgs- W-W and Higgs- Z-Z vertices are both 
proportional to the coupling gHiW, defined as follows [9]: 

9HiW = c pOu + sp0 2i • (A. 3) 

Finally, expressions for the trilinear Higgs couplings gHiHjH k are given in Ref. [9] 
and references therein. 

We now turn to a consideration of the full expressions for the amplitudes in 
Eqs. (3)-(6). The parameters Bf x f 2 always include the Higgs-mediated pieces, but 
also include other contributions with similar kinematical structures. In the text 
these have been parameterized as 

*hh = BfT + 5B hh . 

The expressions for the Higgs-mediated pieces, Bj^ s , are given in Eqs. (8)-(ll). 
In the expressions below we employ the Breit-Wigner form of the propagator for 
sbottom- and stop-mediated graphs, defining 

(V,m 2 ,r) = ^— — , (A.4) 

where we have used a tilde to distinguish the Breit-Wigner propagator from the 3x3 
Higgs propagator matrix defined in Eq. (12). We also use a Breit-Wigner-type prop- 
agator for Z-mediated graphs. Working in Unitary gauge, we take the propagator 
to be %D (p^,m 2 z ,Tz) (—g al3 +PzPz/ m z)- We employ the following definitions for 
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the various combinations of momenta that appear in the propagators, 

M 2 = ( Pl +p 2 ) 2 , 
P 2 = (pi+Pt,) 

£ 2 = (p2+p- tl )' 
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where we note that £ 2 can be written as a function of M 2 , p 2 and various masses. 

The Feynman rules used to build many of the following amplitudes were extracted 
from Ref. [15] (with appropriate modifications to allow for scalar-pseudoscalar mix- 
ing among the Higgs bosons). Amplitudes corresponding to the charge conjugated 
decays can be obtained from those listed below [and in Eqs. (8)-(ll)] by taking the 
complex conjugate everywhere except in the propagator functions Dij and D. This 
is equivalent to making the replacements gHit*t k ~ * 9h t*t an< ^ ^* m ^ ne 

^ j 

various expressions, as well as i — > — i in Eq. (11). 



1. t 2 — ► t 1 H 1 Hi 

The amplitude for t 2 — > t\HiHi is given by 



1 - - sin 2 9 



w 



(O n ) 2 - (0 21 ) 2 + (0 31 ) 2 (4 - c 2 )] t]hh , (A.5) 



where 



__3__Tjt*TTt 

4 cos tV 



12^11 • 



(A.6) 



2. t; 



t^zz 



The amplitudes for t 2 — > are as follows: 

16 



Czz 



Bfz SSS + (2 - y sin 2 ^) ^ , 



- 4 E 1^: 



= -4^(1^ 



2 4 
" 3 

2 4 
~ 3 



sin 2 ^jL>(e 2 ,mf.,r f .)77 ZZ , 
sin 2 6 W ) D(p 2 , m\. , r t ~.) r? zz , 



-^zz — Czz + -Dzz • 



(A.7) 
(A.8) 

(A.9) 
(A.10) 



Each of the non-Higgs-mediated pieces is proportional to the same complex quantity 
Vzz- 



Vzz = Vhh = 



i 77**77* 

4cos 2 ^ 12 11 



(A.ll) 
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3. iz -> i^W+W- 



The amplitudes for t 2 — > t{W + W are given by, 



5 



5 



Higgs 



l+(p 2 -e 2 )^( M2 ' m ^^ 



Cww = 4£> (M 2 ,m 2 z ,Tz) rjww , 



WW 



— 4 



l> (m 2 , m|, r z ) + E |4 1 2 \V tb \ 2 D (p 2 , m?. , r 6 . 



(A.12) 
(A.13) 

rj W w ,(A.14) 



£ w = C w + D w = -4X;|4| 2 |K 6 | 2 ^(p 2 ,m?.,r s .) Ww . (A.15) 

In the numerical work we set V t b = 1. Note that each of the non-Higgs terms is 
proportional to the same complex quantity rjww, which is defined as follows: 



Vww - -^U%U\ X . 



(A.16) 



The amplitudes for £ 2 t 1 ZH 1 are as follows: 

= ~ J; « - ml - ml) D (M 2 , m|, r z ) tol w to 



Bzh 



vg 



m z 

E (^4 - ^sin 2 ^) D (e 2 ,m? ,!>,) g Hli;il , (A.17) 



cos 6> W _ 

C Zi j = 4m z D (M 2 ,m|,T z ) fi^yy to 

E - 1 - 2 <M 2 ') 5 (f 2 . ™|. r £) ) 9Hl£ .«, 

rt* 



" |sin 2 ^ y ) ^ (P 2 ,-I,,r,j g Hliiij , (A.18) 



where 



o 2 - - 
to = Vhh = j \ Q U&Uti ■ 
4 cos^ Vw 



(A.19) 
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